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' Abstract. It is argued that jet modelers have given insufhcient study to the 

^ \ natural magneto-static configurations of field wound up in the presence of a 

■ confining general pressure. Such fields form towers whose height grows with 

' each twist at a velocity comparable to the circular velocity of the accretion disk 

\ that turns them. A discussion of the generation of such towers is preceded by a 

i brief history of the idea that quasars, active galaxies, and galactic nuclei contain 

^SJ ' giant black holes with accretion disks. 

T— I ■ 
> ■ 

o : 

(>0 , 1. Introduction and History 

^ : 

I Before a child learns to run, it must learn to walk. Before it learns to walk, 

' it must learn to stand. I believe that those of us studying jets and winds from 

Q \ rotating objects have been guilty of attempting to run before we have understood 

how to stand. 

Here I give an account of the natural sequence of static magnetic structures 
that are obtained by the continued twisting of a local poloidal magnetic field 
2 , by an accretion disk. My contention is that had we understood such structures 

' long ago, then we would have considered the jets seen in nature to be an obvious 

, consequence of the twisting of the field by an accretion disk in the presence of an 

' external pressure. If we imagine the accretion disk as turning all the time, then 

• i-H . we have to consider a sequence of equilibrium models in which the total twist 

^ ' angles increase linearly with time. This sequence of models gives towers whose 

. height increases with time at a velocity of the order of the maximum rotation 

" " ' velocity in the disk. 

Although I gave the first account of these growing towers in Lynden-Bell 
(1996), the magnetohydrodynamic wind modelers were already much more so- 
phisticated and it is only recently that Li et al. (2001) decided that such an 
elementary model may give a fundamental clue to what is really going on. Their 
numerical method was not capable of exploring such models beyond the first 
two twists. While these are highly indicative of what happens, the many-twist 
picture is best understood from the somewhat more realistic version of my 1996 
paper given here. The essentials of the process can be understood without any 
detailed mathematics and with only the minimal knowledge of magnetohydro- 
dynamics contained in Alfven's theorem that magnetic flux moves with the fluid. 
Those primarily interested in the magnetic towers model of jets may be advised 
to move direct to §2 unless they have an interest in the history of the subject. 

I was originally invited to give an account of black holes in galactic nuclei, 
a subject I had some hand in starting, but by now there are those far better 
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qualified to discuss the recent exciting data and David Merritt has given a mag- 
nificent critical account here. My counter-suggestion of giving some account of 
the history of black holes in galactic nuclei followed by my recent work on jets 
was accepted, so I shall now turn to history. 

Since the Isaac Newton Group is running this conference there is no better 
starting point than Newton's first Query in his Opticks of 1704: 

And do not Bodies act upon Light at a distance and, by their action, bend 
its Rays, and is not this action (caeteris paribus) strongest at the least 
distance? 

In 1784, moving on by 80 years we get to the Reverend John Michell's great 
paper that not only predicted giant black holes but told us how massive they 
would be and how they would be found. Michell was a remarkable scientist, 
one of the founders of the geology of the UK through his early stratigraphy. He 
was the first to suggest that earthquakes traveled as waves. Other achievements 
of his were the experimental demonstration of an inverse square law between 
magnetic poles, the determination of the first luminosity function (using the 
Pleiades), the statistical demonstration that binary stars were gravitationally 
associated, a method for measuring the motion of the Sun through nearby stars 
and the invention of the apparatus for measuring G, the gravitational constant. 
This was perfected and first successfully used after his death by his great friend 
Henry Cavendish. Erasmus Darwin, grandfather of Charles, was among Michell's 
more prominent students at Cambridge. Michell in a paper in the Philosophical 
Transactions of the Royal Society for 1784 argued that bodies of greater than 
500 times the diameter of the Sun and not of lesser density would so attract the 
light that it could not escape, rendering such great masses invisible to our senses. 
Nevertheless, he thought, they might be detected through observations of small 
satellites circling about them. It is interesting that (500)^ Mq = 1.25 x 10^ Mq 
(Mitchell 1784), while the first definitive observation of a giant black hole by 
Miyoshi et al. (1995) used material circling around it to determine its mass — 
3.6 X 10'' Mq. Such predictions more than 200 years in advance of their time 
are extremely rare! 

We know that Michell's idea traveled to the French academy not only by 
the customary exchange of journals but also it was discussed in an exchange 
of letters between the botanist Sir Joseph Banks, President of the Royal So- 
ciety, and Benjamin Franklin, who was then (1783) US Ambassador in Paris. 
Two years after Michell died, Laplace described the idea of dark stars (without 
attribution) in his Systeme du Mond^ Laplace (1795), using exactly Michell's 
argument translated into French, but in 1802 Young published his two-slit ex- 
periment demonstrating the wave nature of light (Young 1802). This put in 
doubt Michell's arguments, based on the corpuscular theory of light, so Laplace 
deleted the passage on dark stars from the second edition of his book. 

The well known historian of science Agnes Clarke, writing an account of 
Michell's life and work for The Dictionary of National Biography (Clarke 1917), 
ended her list of Michell's scientific achievements with the words 



^Available in English translation (Laplace 1809). 
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But [he] speculated fruitlessly on a supposed retardation of light through 
the attraction of its corpuscles by the emitting masses. 

This almost certainly reflects the received wisdom of the British astronomers of 
the time. 

Returning to the 1780s, Sir WiUiam Herschel (1786) remarked upon the 
bright nuclei of the galaxies we now know as NGC 4151 and NGC 1068 and 
in 1850 Lord Rosse, using the giant telescope he had developed at Birr Castle, 
discovered spiral structure in M51 and other galaxies (Rosse 1850). 

Coming to the 20th century, Vcsto Sliphcr (1917), taking long time-exposed 
spectra photographically with a 36 inch telescope, discovered broad emission 
lines in these galaxies. They were so broad that he decided some non-Doppler 
broadening mechanism must be involved. It is seldom pointed out that most of 
the early redshifts of galaxies are due to Slipher's work. He did not claim that 
the Universe was expanding because he thought similar work in the southern 
hemisphere might find almost nothing but blueshifts. However, his work was 
used by Eddington (1922) in his book on general relativity. 

In 1918 Curtis photographed the jet at the centre of M87 (Curtis 1918). 
Even the emission mechanism remained a mystery until Baade showed it to be 
polarized synchrotron emission in 1952. Synchrotron emission was first worked 
out by Schott (1912) but it was not recognized as an emission mechanism in 
astronomy until the 1950s. 

The major postwar development of astronomy came out of British develop- 
ments in RADAR. In the years 1955-1975 Sir Martin Ryle led his Cambridge 
team to a series of remarkable discoveries which impacted almost all parts of the 
subject from stellar death to cosmology (Ryle 1968). Thus the third Cambridge 
Catalogue of radio sources and their 3C numbers became deeply embedded in 
the subject. Ryle and Hewish together developed aperture synthesis as a tool. 
In this they were helped not a little by a secret weapon. David Wheeler in 
Cambridge had developed a program for inverting Fourier transforms which was 
many times faster than others then in use. When some ten years later the 
method was rediscovered and published by Cooley and Tukey it rightly became 
famous. 

Henry Palmer of Jodrell Bank led a team that developed linked radio tele- 
scopes over considerable distances and showed that a small fraction of the 3C 
sources contained small-angular-diametcr sources down to a few arcseconds. 
Hazard realized that an accurate position for one of these could be obtained 
from a lunar occultation that would take place in Australia (Hazard, Mackay, &: 
Shimmins 1963). Such positions were vital for the identification programs which 
Schmidt and Sandage had started at the Mt Wilson and Palomar Observatories. 
The first quasar spectrum to be taken was that of 3C48 by Greenstein. This 
showed emission lines at wavelengths that defied identification. Later that year 
(1963) Hazard sent Schmidt the lunar occultation position of 3C273. Luckily 
this showed two lines whose wavelength ratio exactly matched those of the hy- 
drogen spectrum at the amazing redshift of 2; = 0.158, quite unheard of for a 
12**^ magnitude object. 

Schmidt's (1963) identification was soon confirmed by Oke's (1963) detec- 
tion of a third hydrogen line in the infrared. With the clues provided by 3C273, 
Greenstein's spectrum of 3C48 was then interpreted as having a redshift of 
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z = 0.367. Soon afterwards, Greenstein &; Schmidt (1964) collaborated on a 
paper that demonstrated that the high redshifts of the quasars could not be of 
gravitational origin. The emission lines needed a large volume for their forma- 
tion and in a strong gravity field this would involve much greater line broadening 
than the emission lines showed. Hoyle and Fowler tried to get out of this by 
putting the source of emission at the centre of a very massive large stellar system 
which provided the potential but such a contrived model acquired few backers. 

Early attempts to find cluster galaxies associated with quasars failed and 
this let to the idea that quasars were not associated with galaxies at all. Novikov 
(1964) and Ne'eman (1965) independently put forward the idea that they were 
white holes — time-delayed pieces of the Big Bang. It was Salpeter who first 
considered a massive object ploughing through a galaxy and accreting interstellar 
gas. Assuming that the material gradually percolated down to the last stable 
circular orbit, he gave the 5.7% efficiency of turning mass into energy. The 
properties of the orbits in Schwarzschild's metric had indeed been known to 
relativists for some years, Zel'dovich among them. 

Three papers laid the foundation of the thin accretion disk model now 
widely used 

• Salpeter, E.E. 1964, 'Accretion of interstellar matter by massive objects', 

ApJ. 140, 796 

• Lynden-Bell, D. 1969, 'Galactic nuclei as collapsed old quasars', Nat, 223, 
690 

• Bardeen, J. 1970, 'Kerr metric black holes', Nat, 226, 64 

G. Burbidge (1959) gave minimum energy estimates of order 10^^ erg for 
a number of radio sources and rightly emphasized that this was a very large 
number, while well before quasars Ambartsumian had emphasized the extraor- 
dinary behavior and high energy emissions from galactic nuclei (Ambartsumian 

&; Shakbuzian 1958). 

The main points of my 1969 paper were: 

1. 10^"*^ erg weigh y2l0^ Mq. If these ergs arise from nuclear fusion then 
10^ Mq must be involved. But quasars vary in as little as 10 hours and 
G(10^ M0)2/(1O light hour) » 10^^ erg. The assumption that the energy 
is nuclear therefore leads to an even greater gravitational binding energy 
which must have been lost. Thus the assumption is wrong. 

Hence most of the power of quasars is gravitational in origin and masses 
~ 10^ Mq suffice. But there arc no dead states like white dwarfs for such 
masses, except black holes. If emission is not 100% efficient remnants must 
remain. 

2. I next considered the number of dead QSOs — based on Sandage's (1965) 
estimates of them — and obtained 

/ number of clusters \ / number of \ / number of 
Y of galaxies J \ dead quasars J \ galaxies 
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3. How do we hide dead quasars of ~ 10^ Mq when they still gravitate? 
They are likely to be in dense places and to surround themselves with 
stars. Thus the obvious place to find them is in galactic nuclei. 

4. What should dead quasars look like? Taking them to accrete via a friction 
caused by magnetic torques, I derived a temperature distribution T{r) oc 
r~^/^ and adding many such rings of black hole body emission gave me 

Si, oc fI^^v^/^ exp -(/ll^/A^Tmax), 
where Fm is the rest-mass flux down the hole. 

5. Observational predictions were that most large galaxies should have nuclei 
with high M/L when inactive. In calling the paper 'Galactic nuclei as 
collapsed old quasars', I was raising the question "Are galactic nuclei just 
stars gathered around such black hole remnants of quasars?" 

The masses of the black holes estimated in 1969, mainly from Merle Walker's 
observations using the Lallmand image-tube, are compared with modern values 
in Figure 1. It is clear that these old estimates are mostly too high by a factor 
of order ten because no allowance was made for the mass of stars in the nuclei. 
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Figure 1. My 1969 estimates of Black Hole masses o, compared with modern 
determinations - ★ good, and • less definitive. 
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At the time, most working astronomers were highly sceptical of the whole 
picture, but Maarten Schmidt volunteered that it seemed the best theory avail- 
able to explain the amazing phenomena. 

It was clear that as a black hole accreted from the last stable orbit it would 
accrete angular momentum as well as mass. It thus became important to work 
out accretion disks with Kerr holes. I had a friendly competition with Jim 
Bardeen to do this and his knowledge of general relativity gave him a head start 
which he easily kept. His short but beautiful paper (Bardeen 1970) shows how 
a Schwarzschild hole evolves into an extreme Kerr hole as it accretes. Some of 
this is covered in my Vatican review (Lyndcn-Bcll 1971a) and radio observations 
of the Galactic Centre were made then (Ekers &: Lynden-Bell 1971). 

In 1971 I reviewed the data on the Galactic Centre with Martin Rees 
(Lynden-Bell Sz Rees 1971) and some further details of Kerr disks and vortices 
were given later (Lynden-Bell 1978, 1986), but by then the field had become 
overpopulated. 

A maximum likelihood method for determining luminosity functions was 
developed and applied to quasars and miniquasars (Lynden-Bell 1971b, c). 

2. Jets from Accretion Disks 

Jets were not generally acknowledged to feed radio lobes before Rees's (1971) 
theory. Jets are seen in nature associated with radio galaxies (M87, Cygnus A; 
Hargrave Sz Ryle 1974), quasars (3C273, 3C47), young stars (Beckwith et al. 
1990; Bouvier 1990) Hcrbig Haro objects (Bontemps et al. 1996; Burrows et al. 
1996), dying stars (SS 433), and micro-quasars (Mirabel & Rodriguez 1999). 

One common feature is that the jet velocities are of the order of the circular 
velocity that woTild balance gravity at the central object's surface. 

Theoretical models of jets fall into three classes: 

1. Hydrodynamic models collimatcd by Laval nozzles (Blandford & Rees 
1974), or by vortices around black holes (Lynden-Bell 1978), or by self- 
similar thick disks (Gilham 1981; Narayan & Yi 1995; Narayan, Barret, & 
McClintock 1997). 

2. Wind models in some way collimated by the local magnetic field of the 
rotating source (Mestel 1968; Blandford & Znajek 1977; Sakurai 1985, 
1987; Heyvaerts & Norman 1989; Appl & Camenzind 1993; Lery et al. 
1998; Lynden-Bell 1996; Okamoto 1997, 1999; Li et al. 2001). 

3. Models collimated by a large-scale pre-existing magnetic field (Lovelace 
1976; Blandford k Payne 1982; Shibata k Uchida 1985, 1986; Begelman 
& Li 1994; Bell & Lucck 1995; Lucck & Bell 1996; Ouyed & Pudritz 1997; 
Ouyed, Pudritz, & Stone 1997; Vlahakis & Tsinganos 1998). 

Models in class 3 certainly work but their collimation is no surprise since 
the large-scale magnetic field provides collimation at great distances. A body 
of plasma initially fired across the field at some angle 6 to it would have its 
transverse motion resisted by the field (a high conducting background medium 
being assumed) but its cos 9 motion along the field would continue. Eventually 
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it would be collimated to move along the field. Thus, to some degree, all such 
models provide collimation by fiat in the boundary conditions. They cannot 
provide collimation in dying stars or bodies such as SS 433 in which the jets 
precess (Margon 1984). 

The prevalence of electromagnetic phenomena associated with all jets and 
the acceleration of some particles to very high energies with relativistic motions 
strengthen the idea that the jets themselves may be an electromagnetic phe- 
nomenon and even the early advocates of purely hydrodynamic jets have more 
recently turned to magnetohydrodynamic models. This is partly because of the 
difficulty of getting sufficiently high speeds but also because it was difficult to 
obtain the extremely narrow jets like those seen in Cygnus A and HH 30. 

Here we shall therefore concentrate on the models in class 2, in which the 
jets are collimated not by external fields but by fields associated with the objects 
themselves. Work in this area began from consideration of stellar winds and the 
angular momentum that they transported thereby braking stellar rotation, but 
it took on a new life in the discussion of pulsars. Following Sakurai's early work 
there has been much discussion of asymptotic collimation of such winds, however 
the whole subject has been put in a well ordered form in Mestel's (1999) book 
so the interested reader may find it there. I gave a brief review of the still active 
controversies in my introduction to the Royal Society's meeting on Magnetic 
Activity in Stars, Disks and Quasars (Lynden-Bell et al. 2000). 

As stated in the introduction, we shall now turn to my growing-towers pic- 
ture (Lynden-Bell 1996) based on motion through a sequence of static structures 
that inevitably occur when a poloidal magnetic field is progressively twisted. 

2.1. Magneto-Static Theorems and Rough Estimates 

We consider a magnetic field anchored on the plane z = and confined by an 
ambient external pressure po due to an ionized medium. Where there is field 
there is no gas pressure and except at the surface of the volume V occupied by 
field (and z = 0) the field configuration is force free. The total energy of the 
configuration may be written 

W = Wr + W^ + W, + Wp = ^J [Bl + Bl + BfjdV + poV. (1) 

We have taken cylindrical polar components of the field for our later convenience 
but the theorems proved are not confined to fields of any particular symmetry. 
Now the field wiggles its way to the minimum energy configuration subject to 
the constraints imposed by the boundary conditions on the disk. These define 
the vertical field component B^ on z = and the twist angles of each field line 
which are most easily conceived for the axially symmetric case in which they are 
measured as twists about the axis. It may be shown that the force free equations 
follow from minimizing the energy subject to the constraints. 

Theorem 1 For any finite magneto-static configuration anchored onz = then 
on any cut z = constant 

WR + W^ = Wz+ Wp, (2) 

where wr ,^ ^ = [1 / (Stt)] J B^ ^ ^dA, Wp = poA and A is the area of the volume 
V intersected by the cut. 
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Figure 2. A slice of a force-free field structure. The twists about the axis 
are not shown. The theorem is true without assumptions of axial symmetry. 



Notice that the capital W are related to the small w by integration over z. 

Proof Consider the slice between the cut at z and one at z + dz. Lift the whole 
field configuration above z + dz rigidly by an amount (A — l)dz, expand the slice 
to thickness Xdz and leave the configuration below z unchanged (see Figure 2). 
In the vertical expansion of the slice flux conservation gives 



Br Br/ a , B(j, ^ B^l A , B^ 
Thus in terms of the old fields and old w 



, dAdz XdAdz. 



AW 



(WR + w^) (A ^ - 1) + (w, + Wp) (A - 1) 



dz. 



Now the original situation must be a minimum energy one, so d{AW)/dX must 
be zero at A = 1. This condition yields at once wr + = Wz + Wp (QED). 

Theorem 2 For any finite magneto- static configuration anchored on z = 

Wz = Wp+y2Wo, (3) 

where Wq = [1 / (-^^)] / BzBpjRdA evaluated on z = 0. 

Proof Consider a lateral expansion of the volume V occupied by field such 
that R — > fj,R, z ^ z. Conservation of flux across elementary areas Rd(l)dz, etc., 
gives Br Br/^i, B^/fi, B^ B^/fi^, dA n^dA 

AW = Wz{fi-^-l) + Wp{fi^-l). 

Now we must take care because, unlike our former transformation, this one moves 
the anchored foot-points on z = 0. Luckily the principle of virtual work comes 
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to our aid. We have to add to AW the work done in infinitesimal movement 
S/iR of the foot points by the forces of constraint, viz. BjiBzdA/{4:Tr), so 



dAW 



+ / BnB.RdA = 0. 



Thus we obtain 



W, = Wp+ V2W0 



(QED). To those who beheve that the pinch effect needs no external pressure 
it comes as a surprise that SttVF^ = / B^ji^^fi^dAdz, which is independent of 
fj,. Thus there is no tendency for overall contractions in a purely toroidal field. 
It is perhaps of interest to minimize the energy of a toroidal flux F contained 
between cylinders of radii Ri and R2 of height Z. One finds B^p on 1/R: 

1 F^Z-^ 



4\n{R2/Ri 



Thus if i?2 and Ri are increased by the same factor there is no change in VF^, 
in agreement with what we found above, but if R2 is fixed then decreases as 
Ri is reduced giving rise to a pinch. In fact B'^oaX/R^ shows that the toroidal 

magnetic field acts as a pressure amplifier delivering on Ri a pressure (i?2/^i)^ 
times that exerted on R2. The whole pinch effect is a pressure amplifier — without 
anything to push on outside it has nothing to amplify and indeed the field would 
expand outwards as well as inwards. 

To orient ideas it is useful to make some very crude dimensional estimates. 
Suppose the typical radius out to which the field expands laterally is R2 and we 
deal with an axially symmetrical field configuration which is wound up by the 
turning of the accretion disk. Let the field pass upward through the disk at small 
radii and return downward at larger radii and let the total upward poloidal fiux 
be F. With each turn of the upward fiux relative to the downward a toroidal 
flux equal to R will be generated, so after n turns the toroidal flux will be nF. 
If Z is the total height of the configuration then the volume V occupied by field 
will heV = -kR^Z a typical value of B^ will be nF/ {R2Z), a typical value of Br 
will be F/{^/27^RZ), and a typical value of B^ will be 2F/{ttRI); the ^/2 in Br 
comes from the fact that the typical nR for a radial field line is 27r(i?2/2)^/^, 
while the factor 2 in the B^ comes from the fact that the flux goes both up and 
down. 

Dividing (2) by A/ {Sir) yields the interesting exact result 

{bI) + {bI) = {b^,) + 87tpo, (4) 

where the averages are taken over any z = constant plane. Putting in our rough 

estimates we get 
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Now consider the case of no bounding pressure. Then 

i?2 



| = i(i + 2.V), (0) 



SO Z/R ^ (7r/2)n for n large, which appears to show that the colhmation Z/R 
increases with every turn n. On first finding this I was much encouraged and 
immediately turned to get a more exact model. When I solved it 1 was deeply 
disappointed because before the field has twisted by even one turn it relieves 
itself by expanding to infinity with an opening angle of 120° — very far from the 
sort of colhmation we seek. Figure 3 shows the field configurations for twist 
angles of 130° and 180°. When 208° = 27r/\/3'' is reached the field is outwards 
along cones up to 60° to the equator. This disconnection at infinity is precisely 
what was predicted by some pretty theorems of J.J. Aly (1995) but, although it 
provides a beautiful illustration of these, it fails to give the accurate colhmation 
we seek. However, those interested in galaxies with conical outflows should 
perhaps pay more attention to this magneto-static result described in detail in 
Lynden-Bell & Boily (1994). I was so disappointed in the result, which was 
initially found in 1979, that I laid the whole subject aside and only returned to 
it in 1994. It took two more years before I realised that victory might still be 
snatched from the jaws of defeat. Expression (6) failed to give colhmation only 
because n could not be made large without the field disconnecting at infinity. If 
a way could be found to prevent the field reaching out to infinity then perhaps 
n could be made large after all. This led me to introduce the ambient pressure 
Po to confine the field. It then takes an infinite amount of work for the field to 
extend to infinity so this will not occur at finite twist angles. Let us now turn 
to equation (3) and insert our rough estimates. On dividing by V/Stt we find 

87rpo + 2V2— 2-2, (7) 



^2^4 ""^u ttRIz^' 

where we have left the Bz field on the disk at B^o since this is fixed but used 
our estimate of Br there too. Eliminating po between (7) and (5) we find 



|; = ^(2V2 + 1 + 2.V), 



so asymptotically 



just a factor ^/2 less than our first wrong estimate without pq. With Z^/R'^ 
large the last term in (7) is negligible compared with the first, so Wq may be 
neglected compared with Wz- Furthermore, with n large Br becomes negligible 
compared with B^f,. Thus highly wound configurations become tall towers. For 
such we have from (4) at each height 

Bl) = (bI) + Sttpo. (8) 

Furthermore, we may consider local transformations that expand R ^{z)R 
where ii{z) is slowly varying but becomes one except in a region around some 
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chosen height. Even after this transformation the fields will still have Br small 
so its square is negligible compared with B^. Prom such transformations, which 
leave Wq unchanged, we can deduce that at any great height (3) may be replaced 
by Wz = Wp, so on dividing hy A/Sw 

(bI) ^ Svrpo. (9) 

Hence from (8) we have everywhere well above the base z = 

y2{Bl) = {Bl) = 8Trpo, (10) 

where the averages are taken over areas at constant z. 



3. Better Estimates of the Field Structure 

Consider the tube of force that intersects the accretion disk in the circle of 
radius Ri. Let the poloidal flux rising within this tube be P. Then we may 
label each field line by the value of P and indeed Bz = +[l/{27rR)]dP/dR, 
Br = —[l/[27rR)]dP/dz. If this field line returns to the disk at some larger 
radius R2 then there will be a differential twisting due to the fact that on the 
accretion disk fli > ^2- We define ^{P) = i^i — Q2 so ^{P) is the rate of 
twisting of the field line labeled P. Now suppose that at any given time t when 
the twist angle will have accumulated to be $(P) = Cl{P)t the field line labeled 
P rises to a maximum height Z{P). Then the total twist per unit height will 
be $(P)/Z(P). Now one turn (a twist of 27r) will cause a toroidal fiux between 
P and P + dP of just dP so the toroidal flux per unit height between P and 
P + dP will be 

27rZ(P) ^ ^ 

Now only those fleld lines with P' < P will reach to height Z{P) so the total 
toroidal flux per unit height at Z will be 



B^dR= / „ l,L dP' 



k ^ Jo 27rZ{P') 
so 

where P is the flux label of the line that just reaches height Z{P) and no further. 

The radius of the area occupied by magnetic field at height Z we call R{Z{P)), 
or R{P) for short. This is not the radius at which the line labeled P reaches 
height Z. 

Now the flux P is all the poloidal flux that crosses height Z once on the 
way up and once on the way down, so the average {\Bz\) is 2P/ ^'K[R{Py<^^. 
bmce Bz = {2TrR)-^dP/dR we may write 

, , rPiP) / r ^\2 

{bI) = (27ri?)-i {dP/dRf dR = (2P/ {vr [R{P)f]) I, (13) 
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Figure 3. Field lines splay out and reconnect at infinity after a twist of only 
27r/\/3 rad sa 208° unless confined by an ambient pressure. 



where I is the dimensionless integral given by 



~\ 2 ~ 



and 



I= - I (dP/dR) R-'dR 



P = P'{R')/P, R = R'/R{P). 



(14) 



(15) 



By the theorem that the mean square is greater than or equal to the square of 
the mean modules / > 1. 

Now in a self-similar distribution of field / would be independent of Z so 
we shall suppose that at heights well above the base at z = 0, / depends on Z 
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only weakly if at all. Using (13) for (^B^j in (10) we find 

n[RiP)f = P[I/i27rpo)]'/\ (16) 
so neglecting the variation of log /, 

dlnii:(P)/dlnP= 1/2. (17) 

Now let (bI)/bI = J. Then from (10) and (12) 



/„ 



Differentiating with respect to P and using (17) and (16), 

$(P)/Z(P) = 47r{7rpo/J)'/'R{P)/P = {^I / jf'^ /R{P), 

where we have neglected the variation of J with height. Hence 

Z{P) _ $(P) 
R{P) ~ {81/ jy/^' 

or, since #(P) = 0(P)i, 

Z{P) = {8I/jy/^n{P)R{P)t, (18) 

so the maximum height of each field line grows linearly in time at a velocity 
related to the disk's circular velocity, and the area occupied by the field at 
height Z{P) is proportional to P as given by (16). 

To proceed further we need to determine 0(P), which is one of the inputs 
of our model since it is determined by the distribution of the poloidal flux over 
the disk as well as the disk's rotation. 

We shall take our accretion disk to be Keplerian outside some radius Rq so 
that the circular velocity 

V = Vo{Ro/Rif/^ {Ri > Ro) 

however, such a law cannot continue down to Pi = as it would lead to infinite 
velocities. We shall therefore take circular velocities of the form 



V = Vo 



1 + {Ri/Rof 



-1/4 



My original accretion disk model of quasars (Lynden-Bell 1969) relied on shear- 
ing and reconnecting magnetic fields giving stresses to provide the torques which 
allowed the material to lose its angular momentum and be accreted. It gave mag- 
netic fields B oc P~^/^, and more recent accretion disk models (e.g., Shakura 
& Sunyaev 1973, 1976; Tout & Pringle 1992), which rely on the less reliable 
a— viscosity, give the same R dependence for magnetic fields. However, again 
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the field must reach some finite value as Ri becomes small so we shall take on 
the disk z = 

= Boil + {R,/Rof (19) 
This gives a convenient flux function 



P = Fo^{ l + iRi/RoY 



3/8 



for P>Fo 



(20) 



where Fq = ttR^Bq. We shall assume there is a total poloidal flux F so P < F. 
Eliminating Ri/Rq we obtain 

V = Vo[l + (3/8)P/Fo]"^/^ 
or writing J7o = Vq/Po we have 



(for P <F) 



V/Ri = n = no[l + (3/8)P/Fo]"'/' 



1 + 



3P 
8^ 



8/3 ^ -1/2 
- 1 > 



(21) 



where the last factor is Rq/Ri rewritten and we have neglected $1(^2) compared 
with n(Ri). 

For small P,n = no{P/Fo)-'^/^. 

Wc arc now in a position to determine the shape of the magnetic cavity 
occupied in the magnetic field. From (16) P oc [P(P)]^ so we get O in terms of 
R{P) from (18). For R{P) < b 



n = n, 



1 + (3/8) (P(P)/a) 



-2/3 



8/3 



-1/2 



(22) 



where a? = tt ^Fq [//(27rpo)]^^^ , which is constant and = {F/Fo)a?. Thus 
from (18) writing x = R{P)/a we have for x <b/a: 



Z{P) = {8I/jf/'^a^Qt 1 + (3/8)a; 



-2/3 



1 + (3/8)x^ 



8/3 



- 1 



-1/2 



X. 



(23) 



The factors after the first bracket — > 1 for x small. Z[P) gives the height and 
R{P) alias x gives the radius of the cavity occupied by field so a plot of this 
equation gives the shape at each time. R{P) has a maximum at b and at lower 
Z it remains at h down to Z a where the field configuration becomes that of 
Figure 3. The shape of the cavity is plotted as Figure 4. 

The stability of the towers is an interesting subject that cannot be treated 
here. However, some general guidance can be given. Long thin towers in com- 
pression arc unstable to sideways buckling like an Eulcr strut. If by contrast 
they are in tension, they are stable to such buckling. However, even a liquid is 
stable under compression if it is surrounded by other liquid at the same pres- 
sure. Thus if wc imagine a long thin tower of liquid it will only buckle if asked to 
support a pressure in excess of the ambient pressure of the surrounding liquid. 
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Figure 4. The magnetic cavity occupied by the field becomes systematically 
taller with time but its lateral extent remains the same. 



From this example it is reasonable to expect our magnetic towers to buckle only 
if they experience a net compressive force greater than the ambient pressure. 
With ((-B|) — {B1))/S-k — po = our towers are just stable by that criterion. 
However, we have neglected two other effects. Firstly, magnetic field is buoyant, 
so if pq decreases with height this should require a net tension from the magnetic 
stresses. Secondly, ram pressure at the head of any advancing jet will require 
some extra compression along the jet. As these two effects act in opposite di- 
rections we must leave the stability of the towers undecided. Finally, it is not 
clear that all real jets are stable to buckling; that of 3C273 appears bent on 
milliarcsecond scales. 



4. Conclusions 

In reality jets are dynamic whereas the towers we have calculated are static. In 
so far as the ram pressure and other inertial effects can be neglected, evolution 
through growing towers gives a realizable model. But such models which are very 
useful for guidance and understanding, should not be considered as giving the 



Why Do Disks Form Jets? 



227 



exact shapes to be expected from fully dynamic jets. However, they do supply 
an interesting and provocative answer to the question "Why are there jets at 
all?" Theories of jet models that fail to answer that question but assume an 
imposed flux of material from their base may be missing this point. Our picture 
also gives a good explanation of why jets advance at speeds directly related to 
the circular or escape velocity of the inner parts of the disks from which they 
emerge. 
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